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1. Introduction 

In the last few years the study of sohd-state systems of nanonietric dimension in which 
light interacts with matter has become a widely studied subject. The study of excitons 
interacting with a confined mode of light has made possible the observation of two 
different coupling regimes [U E] (strong and weak coupling) . It also made possible the 
observation of collective phenomena of quasiparticles in semiconductor microcavitites 
[31 m [5] . Motivated by very interesting experimental results regarding condensation of 
polaritons [6l [7] it has been proposed that some of the coherence properties of the above 
mentioned system can be understood in terms of an effective pump of polaritons [S', '9] . In 
these references, a finite system of electrons and holes confined in a parabolic Quantum 
Dot (QD) is used to model the matter component of the system. They first obtain the 
dressed states of the finite system model (using a numerical diagonalization procedure), 
in which the Coulomb interaction between charge carriers and the dipole interaction 
between light and matter is explicitly included. Then they use a zero temperature 
quantum master equation in which they include besides the relaxation processes due 
to coherent emission, an incoherent pump of polaritons. The theoretical results they 
obtain, with this model, reproduce the polariton laser threshold reported in pTO] . 
In a recent theoretical and experimental study [11], using a simple Jaynes-Cummings like 
model including dissipative processes, the authors obtain a surprisingly good agreement 
between the calculated and measured polariton spectra. They show that the effective 
dissipative parameters of the system depend on the nominal excitation power density. 
They also discuss the difficulties involved in determining the strong coupling (SC) regime 
mainly because the broadening of the spectral lines hinders the well known anti-crossing 
characteristic feature. The determination of a clear signature of the strong coupling 
regime motivated several recent theoretical works [121 1131 [II]- 

The above paragraph highlights the importance of studying the role of the incoherent 
pumping in the determination of the dynamical regimes in the system. In particular, 
the effect of an incoherent pumping of polaritons or excitons must be clarified. In 
this work we seek to discuss the dynamical effects of the interplay between polaritonic- 
excitonic pumping, by using the simplest model of quantized light matter interaction, the 
Jaynes-Cummings model [151 [IS IE] • We also seek to understand how the entanglement 
between the excitons and photons is affected by the incoherent pumping since we have 
two strongly coupled interacting systems. Some previous works in this direction are 
[18l[l9l[20l[2T]. 

The paper has been written as follows: In section 2 a description of the master 
equation of the system is given. The matrix elements of the density operator related 
to the polariton pumping Lindblad superoperator are obtained. We also show the 
operators and dynamical equations necessary to obtain the emission spectrum using 
the quantum regression theorem. Finally, in this section we derive a closed form 
expression for the Peres positive partial transpose criterion in order to quantify the 
exciton-photon entanglement. In section 3 we compare the effects of both types of 
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pumping in the statistical properties of the steady state density operator. Then we 
show, analytically and numerically, how the polariton pumping affects the emission 
spectra of the system. Finally, in this section it is shown the effects of both types of 
pumping in the entanglement of the system. 

2. Theoretical Background 

We are interested in studying an exciton interacting with the lowest energy (frequency) 
mode of a semiconductor microcavity. The quantum states resultant of the electrostatic 
interaction between holes in the valence band and electrons in the conduction band in 
a solid state system are termed excitons. This quasiparticles might exhibit a discrete 
or continuum spectrum depending on their confinement. In this work we shall consider 
only the lowest energy levels of the system, the ground state, \G) (no excitation, i.e. 
electron in the valence band) and excited \X) state. In this effective formulation all 
the complexities related to the many body problem of considering electrons and holes in 
quantum dot or well are effectively included in the energy separation between the ground 
and exciton states. The light will be treated as a single electromagnetic quantized mode. 
This assumption is subjected to existence of well separated energy modes in the cavity 
we are considering. This assumptions are the ones usually considered in most theoretical 
work [ni[l2l[I71[T8]. 

2.1. Hamiltonian and Dressed States 

The Hamiltonian we shall use to study this system is the well known Jaynes-Cummings 
Hamiltonian |15J: 



H = {(jJx — A)a'a + loxCf cr + 5'(o"a' + a'a), (1) 

where IX),!^) are the exciton and ground states of the matter, a = IG) (X|, cr^ = 
\X) {G\ , 0,0"^ are the annihilation and creation operators of the field, ux is the energy 
required to create an exciton, A = ux — ^0 is the detuning between the exciton and 
photon frequencies, g is the light-matter coupling constant or Rabi constant and we 
have taken h = 1. 

This Hamiltonian is obtained after considering that the exciton is coupled only to one 
mode of the microcavity, that the interaction between them is of dipole type [i.e. 
the spatial variation of the electromagnetic field is small in the spatial dimensions of 
the exciton) and after neglecting counter rotating terms (the so called Rotating Wave 
Approximation) [16]. 

The Hamiltonian H can be diagonalized when written in the basis {\G) , \X)} ® 
{|n)}5!^Q (The Bared Basis). It takes the block-diagonal form: 

{n-l){ux-^)+uJx g\/n \ ,. 

gV^ n{ux - A) i ' ^> 
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when written in the n-th excitation manifold basis {\Xn — 1) , \Gn)}. The eigenvalues 
and eigenvectors can be easily obtained and are given by [16j: 




2nA + A ± V4V + A2) + nujx (3) 



Where A is clock wise rotation matrix by the angle 0„/2 = tan~^ f ^^ j /2. The states 
|n, ±) are the dressed states of the Hamiltonian ([1]). If the states \G) , \X) are the 
excitonic states of the QD then the states \n,±) can be considered the polaritonic 
states of the system. 

The emission spectrum of the system in this idealistic model is given by the tran- 
sitions that can occur between two given dressed states and the values of the transition 
energies are precisely the differences of their respective frequencies times h. 



2.2. Master Equation 

Since the system we are considering is really an open quantum system the effect of the 
environment must be included. The effects of the weak coupling with the environment 
are accounted by writing a master equation in the Born-Markov approximation , for the 
dynamics of the density operator of the reduced (matter-light) system [22]. This master 
equation accounts for the following processes: 

(i) The continuous and incoherent pumping of the exciton. 

(ii) The direct coupling of the exciton to the leaky modes which induces the spontaneous 
emission process. 

(iii) The escape of cavity mode photons out of the microcavity due to incomplete 
reflectance of the mirrors, the so called coherent emission. These photons are the 
ones that are measured to obtain the emission spectrum of the system. 

(iv) An incoherent pumping of polariton (dressed) states. 

The first three processes have been discussed in detail in references [TT|, [HI [T71 [18] 
and we shall omit their details here. An schematic representation of the action of the 
processes involved in the dynamics of the systems is given in panel a) of figure [H The 
fourth process of the above list is intended to cause incoherent transitions among dressed 
states of two consecutive excitation manifolds |n, ±) — )■ \n + 1, ±) as figure ^ panel b) 
shows. This is the equivalent prescription of the polariton pumping considered in [HI |9] , 
in which the dressed states of the light matter Hamiltonian are obtained and are used to 
define raising and lowering operators between two consecutive excitation manifolds of 
polaritonic states. In complete analogy with [H [9] we introduce the following lowering 
operators between polaritonic states: 

P++„ = |n, +) {n + 1, +1 P__„ = |n, -) (n + 1, -I 

P+-n = \n, -) (n + 1, +1 P_+„ = \n, +) (n + 1, -| 
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Figure 1. a) Ladder of bared states for a two level quantum dot coupled to a cavity 
mode. The double headed green arrow depicts the radiation matter coupling g, dashed 
yellow lines the exciton pumping rate P, dotted blue lines the spontaneous emission 
rate 7, dashed dotted red lines the emission of the cavity mode k and grey wavy lines 
the polariton pumping process Pp. b) Transitions due to the polariton pumping term 
in the master equation ([5]) in the ladder of dressed states. 



With this definitions the master equation for the system takes the form: 



— p = i[p,H] + -{2apa) 



dt 



a^ap 



pa^a) 



1 + + 

— {2apa^ — cr^ap 



pa'^a 



(5) 






aa' p — paa 



' r, / J ^^ijnP^ijn ^ijn^ijnP P^ijn^ijn 



n,t,J 



Where k is the decay rate of the cavity photons due to the incomplete reflectance of the 
cavity mirrors, 7 is the exciton decay rate due to spontaneous emission, P is the rate 
at which excitons are being pumped and Pp is the rate at which polaritons are being 
pumped. The indices i and j take the values {+, — } and n E N. The polaritonic bath 
that is associated to the polariton pumping term in the last equation can be thought 
as a resonant coupling between the bare electron states and the intersubband cavity 
polariton excitations |23j . 

For studying the dynamics of the system we shall write the equations of motion of 
the matrix elements taken in the bared basis. The evaluation of this elements is straight- 



forward except for the term Lp^ [p] = | J2n,i,j '^PijnPPi 



ing the fact that |n, +) (n, +| 



n, 



-)(^, 



tjn 



^ijn^ijnP P^ijn^ijn 



\Gn){Gn\ + \Xn-l){Xn 



Us- 
1 1 and that 
the trace of an operator is invariant under unitary transformations it is easily seen that: 

(i, n\ Lp^ \p] I j, m) = SG,iSG,jSm,nPGm-l,Gm-l + k,x5j,x5m,nPGn,Gn (6) 

+ SG,i^G,j^m,nPXm-2,Xm-2 + Si,xSj,xSm,nPXn~l,Xn-l 
~ ^GAPGnjm — Si,XPXn,jm — ^GjPin^Gm — Sj,XPin,Xm, 

which is independent of both g and A. The dynamical equations for the populations and 



coherences in the bare basis are presented in the appendix Appendix A For simplicity 
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in what follows we have taken as unit of frequency (equivalently energy) the Rabi 
constant g = 1 nieV which is a typical value of the light-matter coupling constant 
in semiconductors. All the quantities in what follows are given in units of g. 

2.3. Quantum Regression Theorem and the Emission Spectrum 

One of the few things that can be measured directly from a quantum system is its 
spectrum. To obtain the emission spectrum of a system we need to take the Fourier 
transform of the first order correlation function (a^(t + r)a(t)), which requires the 
knowledge of the expectation value of two operators at different times. To obtain the 
dynamical equation of such correlation function we take advantage of the Quantum 
Regression Theorem (QRT) [21] which states that given a set of operators Oj satisfying, 
f (0,(t + r)) = E, L,, {Ou{t + r)) , then f (0,(t + r)0{t)) = Y.u hk {0,{t + T)0{t)) 
for any operator O. We follow Tejedor and coworkers [TTj and write {a'^it + T)a{t)) = 

J2n ^'^ + 1 ( {'^Gni^ + ''")'^(^)) + {O'Xni'^ + ''")'^('^)) ) ? wherc the following definitions have 
been used: 

al^ = \Gn + 1) {Gn\ (7) 

a^^ = \Xn + l) {Xn\ 

ai = \Xn) {Gn\ 

Cn= \Gn + 1) {Xn - 1\ . 

Note that these operators act between two consecutive excitation manifolds. It turns 
out that the last set operators satisfy the following set of closed differential equations: 

^(4n(^)) = {-P -2P,-i^-nK-^ + lujx) {al^r)) (8) 

+ «:V('^ + l)(n + 2)(a^„^i(r)) 
+ 7(«k(^)) - ^9^/^{Ur)) + igyf^r+l{ai{T)) 

^(4(^)) = «^Vn + l(aJj„(r)) - ig^{a\^_^{r)) 

+ ^(-P - 4Pp - 7 - 2nt, + 2iuJx){(yi{r)) 
+ (n + l)«:((Ti+i(r)) 

^(«k-iM) = ^(4n-i(^)) + «:V^(ri + l)(aL(r)) 

+ ^(-^Pp - 27 - 2zA -2nK + K + 2iujx){a^xn-ii.'^)) 
+ ig\/n + l{Cn{T)) - igy/n{ai{T)) 
(Cn(r)) = -i^v^(aJ.„(r))+i^y^5TI(a5f„-i(r)) 



dr 



+ -(-P - 4Pp - 7 - 42A - 2nK + 2iux){Cn{T)) 



+ ^n{n + 2)K{Cn+i{r)). 
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The QRT implies that the two time operators (ag,^(t + r)a(t)), (a^„(t + r)a(t)), 
(o"^(t + r)a(t)), {(nit + T)a(i)) satisfy equations ([8]), subject to the initial conditions: 

(ajj„(t)a(t)) = VnTTpGn+i,Gn+iit) (9) 

(a^^(t)a(t)) = VnTlpxn+i,Xn+iit) 
{al{t)a{t)) = VnTlpGn+i,Xnit) 
(C„(t)a(t)) = Vnpxn,Gn+i{t). 

The role of the parameters uJx,^,g, K,,'y,P,Pp is twofold, on the one hand they 
determine the dynamics of the two time operators via ([8]), and, on the other hand 
they set the initial conditions ([9]) that shall be propagated according to the dynamical 
equations (IHl). Here we will be interested in the light that the systems emits in the 
stationary limit so that the limit t — ?■ cxo will be taken in equation (Q. 



2.4- Entanglement 

To quantify the entanglement between excitonic and photonic subsystems we use 
the Peres criterion [25]. This criterion says that if the state of bipartite system is 
separable then the eigenvalues of its partial transpose are all positive. For our case 
the density matrices, p, we are considering have only the following non-zero matrix 
elements: pxn,Xn, PGn+i,Gn+i, PGn+i,Xn, Pxn,Gn+i Correspondingly the nonzero elements 
of the partial transpose p^ respect to excitonic subsystem are: 

PGn,Gn = PGn,Gn (10) 

r _ 

PGn+l,Gn+l 
PGn+l,Xn 
PGn,Xn+l = PXn,Gn+l 

The matrix p^ takes a block diagonal form of blocks 1x1 and 2x2, when written 
in the basis {|X0) , |G0) , |X1) , |G1) , \X2) . . . \Gm) , \Xm + 1) , \Gm + 1)} (Notice the 
ordering of the basis) : 

/ Pxo,xo \ 



PGn+l,Gn+l 

r 

PXn+l,Gn 
T 



P 



PGn,Gn PXn,Gn+l 

PGn+l,Xn PXn+l,Xn+l 



\ 



■ (11) 



PGm+l,Gm+l ) 

The eigenvalues of the last matrix are easily obtained. The ones corresponding 
to the left upper and right lower entries (blocks) are pxo,xo and pcm+i.Gm+i and 
are always positive ore zero. The ones corresponding to the 2x2 blocks are 

PGm,Gm + PXm+l,Xm+l ± y {pGm,Gm. — PXm+l,Xm+lY + 4|pXm,Gm+lP j • In Ordcr tO 
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have a negative eigenvalue (and an entangled state) the following condition must be 
met for some n: 

\pXn,Gn+l\ > y/PGn,GnPXn+l,Xn+l- (12) 

Notice that the above inequality can also be obtained by using the criterion recently 
derived in [26] . Then, we can quantify entanglement using the following function, which 
is equivalent to the Peres criterion: 

E{p) = 4y^ (max {O, \pXn,Gn+l\ - y/PGn,GnPXn+l,Xn+l}) ■ (13) 

n 

For Bell-like states pBeii = \ip) {ip\, \ip) = 75 {\Gn + 1) + e*<^ \Xn)), e M, E{p) will 
equal 1. In, particular the polaritonic states \n, ±) in resonance have E{p) = 1. 

3. Results and Discussion 

3.1. Exciton Pumping vs. Polariton Pumping 

In this section, we compare the evolution of some observables in the stationary limit 
as a function of the detuning (A) and the pumping rates of both excitons (P) and po- 
laritons (-Pp). In Figure [21 we present the evolution of the average number of photons, 
(n) = {a^a) in the stationary limit for the case of strong coupling (l/t — 7|/4 ^ (?) as a 
function of A and either P 01 Pp. From this figure several conclusions can be drawn. 
First the pumping of polaritons is more efficient accumulating photons in the cavity, this 
can be understood if one sees the type of transitions that both types of pumpings cause 
in the bare basis ladder of states as it is seen in Figure [TJ On the one hand P causes 
diagonal transitions, whereas Pp induces vertical (and diagonal) transitions that can 
rapidly populate states with a high number of photons. Secondly, and rather surpris- 
ingly, the average number of photons (and as we shall comment later) the observables 
we have monitored are not sensible to the detuning when the term Pp > /t = 0.1. 

We have also calculated the effect of both mechanisms on the population inversion 
{(Tz) and the second order coherence function at zero delay g'^{T = 0) = {o)a)aa) / {o)a) . 
In Figure [2] we present the contour plots of (7^(r = 0) and (o"^). It is clearly seen that 
one can explore different field statistics as P and A are varied and that the population 
inversion grows as a function of P. The corresponding results for the polariton pumping 
mechanism can be summarized as follows: 

• The second order coherence function is almost equal to one, and is independent of 
A for values of Pp larger than k = 0.1. 

• The population inversion presents a similar behavior to 5'^(r = 0). It approaches 
to from below as Pp is increased and it is nearly independent of A for values of 

Pp> K = 0.1 

Comparing the results of both types of pumping it is seen that the polariton pumping 
mechanism is unable to cause a positive population inversion, can greatly increase the 
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Figure 2. Average number of photons, second order coherence function at zero delay 
and population inversion as a function of the detuning A and the exciton pumping rate 
P or the polariton pumping Pp. Parameters: k = 0.1, g — 1, ujx = 1000, 7 = 0, for the 
upper panel Pp = and for the lower panel P = 0. In most experimental situations 7 
is at least two orders of magnitude smaller than the rest of the parameters [27| . For 
this calculations the Fock space was truncated in rimax = 40 



intensity of the light stored in the cavity, and the hght stored in the cavity might have 
Poisson statistics, since the variance of the number of photons equals its mean value 
(although in general it will not be a coherent state since the reduced density matrix of 
the photons will be a mixed state). The first observation is understandable when one 
inspects the action of the term Pp in the bare basis. On the one hand the term Pp causes 
transitions with equal intensity from \X) — > \G) and from \G) — > |X), this explains 
why (az) is close to zero. The fact of being always shghtly negative is related to the 
asymmetry that the vacuum state IGO) introduces between the states G and X, in such 
a way that state vectors with no exciton have a slightly heavier statistical weight. The 
second observation is related to the possibility of having an inversionless polaritonic 
laser ^M\. 



3.2. Emission Spectrum 

Since the effects of the excitonic pumping have been extensively studied in [121 US [E 
in this section we only study the effects of the polaritonic pumping in the emission 
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spectrum of the system. To do so we write the equations of motion of the variables in 
equation ([8]) as: 

-v(t) = A(Pp, P, g, /€, 7, ux, A)v(t) (14) 

where v(t) = {(a^(t)) , (a^oW) , • • • , iUt)) , (aUt)) , (aL-i(t)) , (ak(^)) ,•••}• The 
formal solution of equation ( fT4l) is given by: 

v(t + r) = exp(A(Pp, P, (?, n, 7, w^, A)r)v(t) (15) 

From equations ([8]) one sees that A(Pp, P, ^f, k, 7, wx, A) = — 2PpI +B(P, (7, k, 7, wx, A) 
(I is the identity matrix), i.e. the polariton pumping term is diagonal in equation ([8]), 
and because it commutes with the matrix B(P, g, n, 7, ux, A) it can be factored out in 
equation flTSl) as follows v(t + r) = exp(— 2Ppr) exp(B(P, f;, k, 7, cjx, A)r)v(t) 

The last equation implies that the pumping rate Pp cannot modify the oscillation 
frequencies of the first order correlation function since it acts as a common overall 
decay rate for all the operators involved in equation (jH]). It can only redistribute the 
statistical weights of the different frequencies by modifying the initial values of the two 
time operators that are related to the populations and coherences in equation (J9]). For 
instance in the case where one considers transitions between the ground state \G0) 
and the states \X0) and \G1) only two operators appear in the equations of the QRT, 
{oQoit)) and {ucoit)). This approximation is valid when the pumping (P or Pp) is small 
enough as compared with the losses (7 and k) to not have an average photon number 
of more than one p!8]. In this case the equations of the QRT are: 

(16) 

^ + iujx ig \ f {aloit)) 

-P-2Pp-tA-l+tux J [ (4oW) 

The eigenvalues X± of the square matrix in the last equation will be related to the 
positions u± and widths T± of the emission spectrum {iX± = u± + iT±) and are given 
by: 

X± = -{-3P - 8Pp - -f - 2iA - K + Ucox (17) 

±i^/l6g^^^lP^'^^T2iATKJ^) 

As expected the term Pp only enters as a decay rate that widens both peaks of emission 
in the same way. 

Actually one can go a bit further in the analytical calculation by considering the case 
where there is no exciton pumping, P = 0. In this case one can obtain analytically all 
the eigenvalues of the matrix A = A(Pp, P = 0,g, n, 7, oux, A). To this end notice that 
the only term that couples the group of operators (Cn(^)), {Cn{t)), {O'Xn-ii'^))^ i'^Gni^)) 
and (Cn-i(i)), (o-Lil*^)). («xn-2W)' («Gn-i(^)) ^^ equatlou (^ is the term P {a^cn-i) i^ 

d_ 
(It 




the equation for ^ {cJxn-\li so if P = the structure of the matrix A will consist of 
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block diagonal terms of sizes 2x2 and 4x4 and off diagonal terms over the diagonal of 
the matrix. This implies that to reduce the matrix A to an upper triangular matrix it is 
necessary to rotate each diagonal block and the blocks above it. Once the matrix has an 
upper triangular form their eigenvalues are simply given by the elements of the diagonal. 
Summarizing, the eigenvalues of the whole matrix A are simply the eigenvalues of the 
blocks 2x2 and 4x4. 

The eigenvalues of the 2x2 matrix are given by equation (1T7|) with P = 0. The structure 
of the blocks 4 x 4 is: 

A|4x4= (18) 

' i{uJx — 2A) igVn + 1 —ig^fn ' 

iuJx — '^Pp —igy/n ig\/n + 1 

— I — riK 
ig\/n + 1 —ig^Jn iioJx — A) — 2Pp 

-{n- |)/€-7 



—igy/n igy/n + T i{uJx — ^A) 

y -{n + 1)k - 2Pp 

The eigenvalues of the above matrix are simply given by: 



X±± = -2P„---iA-nK + iux± — ^\/a±Vb (19) 

2 2a/2 ^ 

a= - 8{2n + l)g^ - AA^ + (k - -ff 
b = 256n{n + 1)/ + [{k - 7)^ + AA'^] ^ 

- 16 [{2n + 1)(72 + K^ + 4A2) - 2{2nK + k- 2iA)-f - HAk] g^ 

An interesting question that may be asked is under what conditions the imaginary parts 
of the eigenvalues differ at resonance (A = 0), i.e., under what condition a±Vb < 0|a=o- 
The condition is simply given by: 

, , \l^ — l\ / N 

9 > Va^ 20 

Ay/n 

Which for n = 1 reduces to well the known relation for strong coupling in the absence 
of exciton pumping [16] , but actually also tells that different excitation manifolds might 
be in different regimes so while some excited manifold might be "dressed" while lower 
states might be "bared" [H]. In our case this condition will not be affected by the 
polariton pumping Pp. 

From the above discussion it is clear that the main effect of the term Pp is to cause 
an homogeneous broadening in the emission spectrum of the system. It will also increase 
the intensity of the emitted light since it will increase the population of highly excited 
states in the stationary limit, equivalently the initial values of the two time operators. 
The two mentioned effects can be seen in figure [3] for resonance condition A = For 
low Pp the peaks of the emission appear in u;+ = 1001 and u^ = 999 which are precisely 
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Figure 3. Emission spectrum of tlie system witli different rates of polaritonic pump 
Pp for the parameters tux ~ 1000, k = 0.1, A = 7 = P = 0. For this calculations we 
truncated equations |S] in n„iax — 20 



the transition energies between the states |1, ±) and the vacuum |G0) {ux ± g)- When 
the pumping is increased new hues appear, this hues are associated with the energetic 
transition between the states |2,±) and |1,±) which have frequencies cox ± 5' ± v^g 
(for the parameters used they are approximately Ui ~ {1002.4,1000.4,999.6,997.6}). 
When the pumping is further increased the hues associated with the transition to the 
vacuum become completely shadowed by the widened lines associated with the transition 
|2,±)^|1,±). 



3.3. Entanglement in the Stationary State 



In this section, we use equation ( 11311 to study the entanglement in the stationary 
state. It is important to remember that the Peres criterion is a necessary and sufficient 
condition for having entanglement when the dimensions of the Hilbert spaces (hi) of the 
subsystems considered are dim(/ii) = 2 and dim(/i2) = 2 or 3 [28j, for larger dimensions 
it is only a sufficient condition for having an entangled state. Under certain conditions 
the Fock space of our system can be truncated in one or two photons. This conditions 
are met when we consider small enough pumping rates (P or Pp) as compared with the 
loss rates (k, 7). This hypothesis is numerically con&med by figure H] where it is clearly 
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Figure 4. In the left and central panels we present the fraction of the population 
of the density matrix that has more than two photons as a function of the exciton 
or polariton pumping rates and the detuning A for the parameters k = 0.1, 7 = 0, 
g — 1. The results do not change significatively if 7 = 0.005. The right panel examines 
the entanglement measure E[p) as a function of the exciton pumping rate P and the 
detuning A (The other parameters are the same of the other two panels). It is seen 
that even in the stationary state one can have entangled states, which is not case if 
one only pumps polaritons. The Fock space was truncated in nmax = 40. 



seen that for P or Pp < 0.1 



K + 7 the population corresponding to states with more 
than 2 photons is very small or zero. Then in the blue regions of figure S] we can obtain 
the stationary state of the system considering only a Fock space of 2 photons. The 



exact expressions for a maximum of 1 photon are presented in appendix Appendix B 



One can also find exact expressions for a maximum of 2 photons but they are very 
cumbersome and are no presented here. In the cases where one considers only the term 
Pp (P = 7 = 0) one analytically finds that the state is always separable (truncating in 1 
or 2 photons), this was also numerically confirmed using a larger basis (nmax = 40) for 
the parameter region 10~^ < Pp < 1, < A < 2, k = 0.1 and 7 = {0, 0.005}. When the 
effect of the excitonic pumping is considered we find certain region of parameters where 
there are entangled states (see the right panel of figure Hj). This results highlight the 
very special role of the exciton pumping term (P), because, its variation might induce 
or reduce the degree of entanglement between the exciton and the field mode. 

4. Conclusions 



In this work a mechanism for incoherent pumping of Polaritons was proposed. The form 
of its matrix elements was derived and used in section [2l2] to obtain the dynamical equa- 
tions necessary to propagate the density operator of the system and obtain its emission 
spectrum. The effects of the new mechanism were compared with the effects of the 
exciton pumping. The physical origin of the new mechanism is still controversial and it 
is necessary further work in this direction. In the literature it has been considered an 
effective pump of excitons and photons. The fit that has been done in reference [11] has 
shown that indeed both terms are able to account for a significant amount of the physics 
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in such systems. Nevertheless using out off resonance excitation in quantum wells (QW) 
the group of J. Bloch obtains a polariton laser ^U\, that has been successfully modeled in 
reference j,9j including a effective polariton pumping term, that is essentially equivalent 
to our model. Additionally in reference [6] it has been shown that for QW an effective 
resonant pump to the lower polariton branch is a condition for thermalization of a BEC 
of polaritons. 

In this work it was shown that the polariton pumping term is not able to cause 
a positive population inversion and that above certain threshold drives the field to a 
Poisson-like statistics where g'^{T = 0) = 1. It was shown that the term Pp is unable 
to change the dynamical regimes of the system, and that its effect on the emission 
spectrum is twofold: it causes a homogeneous broadening of the peaks and an increase 
in the integrated emission. Also in this work, it was examined how the incoherent 
pumping mechanisms affects the entanglement of the exciton and the cavity in the 
stationary state. It was shown that the exciton pumping term does not completely 
destroy entanglement as the polariton pumping does. Finally, in section 12.41 a useful 
rule-of-thumb based on the Peres criterion was obtained to determine if the state of the 
exciton and the photons is entangled. 
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Appendix A. Dynamics of the density matrix elements 

The dynamical equations for the populations and coherences in the bare basis are given 
by: 

-T:PGn,Gn = " PpGn,Gn + K, {{n + l)pGn+l,Gn+l - npcnfin) (A.l) 



dt 



+ Pp {PGn-l,Gn-l — '^PGn,Gn + PXn-2,Xn-2) 
+ wVn {pGn,Xn-l - PXn-l,Gn) + lPXn,Xn 



-7:PXn-l,Xn-l — P pGn-l,Gn-l + /« {npXn,Xn — {u — l)pxn-l,Xn-l] 
+ Pp {PGn-l,Gn-l + PXn-2,Xn-2 — 2pXn-l,Xn-l) 

+ 'igyn{pxn~l,Gn — PGn,Xn-l) — lPXn-l,Xn-l 

d 

-T:PGn,Xn~l 

dt 



igy/n{pGn,Gn - PXn-l,Xn-l) + K^/n^U + l)pGn+l,Xr 



2P„ - - + iA - riK + - 

2 "^ 2 2 



+ ( -— - 2Pp - - + iA - UK + - ] pGn,Xn-l 
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Appendix B. Stationary state truncating in one photon 

The populations and coherences of the stationary density operator, truncating the Fock 
space in one photon, are given by: 

Pxo,Gi = i2tg{P + 4Pp + 7 - 22A + k){{P + Pp)K - P^^)) /N (B.l) 

Pgo,go = (4(P + 2Pp)(7 + «:)(P + 4Pp + 7 + K)g' 

+ Pp-f{P + 2Pp + 7 + /€) (4A2 + (p + 4Pp + 7 + Kf))/C 
Pxo,xo = (4(P + 2Pp)(7 + ft:)(P + 4Pp + 7 + K)g^ 

+ (P + Pp)K{P + 2Pp + 7 + /€) (4A2 + (p + 4Pp + 7 + Kf))/C 
PGi,Gi = (4(P + 2Pp)(7 + «:)(P + 4Pp + 7 + K)g^ 

+ Pp7(P + 2Pp + 7 + /€) (4A2 + (p + 4Pp + 7 + Kf))/C 
Pxi,xi = {^9\P + 4Pp + 7 + /s:)(P + 2Ppf 

+ Pp{P + Pp){P + 2Pp + 7 + /€) (4A2 + (p + 4Pp + 7 + fi:)2))/C 
iV =4(P + 2Pp + 7 + fi:)(P + 4Pp + 7 + /s:)/ 

+ {P + Pp + 7)(Pp + /t) (4A2 + (p + 4Pp + 7 + Kf) 
C =N{P + 2Pp + -f + K) 
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